global xx yy
dir1 = ('./figure/');
dir2 = ('./figurew/');
for i = 1:nts
    if restart == 1
        i = nts_start+i;
    end
    load([saldir 'dat',sprintf('%4.4d',i)]);

%     colormap jet
%     pcolor(yy(2:end-1),xx,squeeze(wi(:,2:end-1,1))); colorbar('Fontsize',15); shading interp;
%     title(['vorticity,t=' num2str(i)],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
%     drawnow
%     print(gcf,'-dpng',[dir1,sprintf('w%d',i)]);
%     close
%     
%     colormap jet
%     pcolor(yy(2:end-1),xx,squeeze(phii(:,2:end-1,1))); colorbar('Fontsize',15); shading interp;
%     title(['phi,t=' num2str(i)],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
%     drawnow
%     print(gcf,'-dpng',[dir1,sprintf('phi%d',i)]);
%     close
    
    colormap jet
    set(gcf,'position',[500,500,1000,200])
    pcolor(yy(2:end-1),xx(1:end),squeeze(T0i(1:end,2:end-1)+Ti(1:end,2:end-1)-1)); colorbar('Fontsize',15); shading interp;
    title(['T,t=' num2str(i)],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
%     caxis([-2 2])
    drawnow
    print(gcf,'-dpng',[dir1,sprintf('TT%d',i)]);
    close
    
%     colormap jet
% %     set(gcf,'position',[500,500,2000,400])
%     pcolor(yy(2:end-1),xx(1:end),squeeze(Ti(1:end,2:end-1))); colorbar('Fontsize',15); shading interp;
%     title(['T,t=' num2str(i)],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
% %     caxis([-2 2])
%     drawnow
%     print(gcf,'-dpng',[dir1,sprintf('TT%d',i)]);
%     close
    
%     colormap jet
%     pcolor(yy(2:end-1),xx,squeeze(T0i(:,2:end-1))); colorbar('Fontsize',15); shading interp;
%     title(['T0,t=' num2str(i)],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
%     drawnow
%     print(gcf,'-dpng',[dir1,sprintf('T0%d',i)]);
%     close
    
%     fig = figure;
%     plot(xx,Ti(:,40));
%     title(['T,t=' num2str(i)],'FontSize',15); xlabel('x','FontSize',15);
%     frame = getframe(fig);
%     img = frame2im(frame);
%     imwrite(img,[dir1 'plot' num2str(i) '.png']);
%     close
       
%     colormap jet
%     pcolor(yy,xx,squeeze(T0i(:,:,1))); colorbar('Fontsize',15); shading interp;
%     title(['T0,t=' num2str(i)],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
%     drawnow
%     print(gcf,'-dpng',[dir1,sprintf('T0%d',i)]);
%     close

%     colormap jet
%     pcolor(yy,xx,squeeze(NL(:,:,1))); colorbar('Fontsize',15); shading interp;
%     title(['NL,t=' num2str(i)],'FontSize',15); xlabel('y','FontSize',15); ylabel('x','FontSize',15);set(gca,'FontSize',15)
%     drawnow
%     print(gcf,'-dpng',[dir1,sprintf('NL%d',i)]);
%     close
end

% for j = (nts-50):nts
%     load([saldir 'dat',sprintf('%4.4d',j)]);
%     rms = sqrt(mean(Ti.*Ti,2));
%     plot(xx,rms);
%     
% %     plot(xx(1:(nx+1)/2),rms(1:(nx+1)/2));
%     hold on
% end
% print(gcf,'-dpng',[dir1,'rms']);
